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1. INTRODUCTION 

Nowadays, harvesting energy (EH) from green environmental solar, wind, geothermal, and 
mechanical, the radio frequency (RF) signals can be considered as the prospective energy source in the 
future. RF signals can be proposed as the role of information transmission or harvesting energy in the 
communication network [1-8]. In a traditional network, the security problem to protect the wireless 
communication network is based on cryptographic technologies. Nowadays, the physical layer security (PLS) 
with the new technology is considered as an extensive research direction to protect the modern 
communication network. PLS technology is considered as a novel way to protect the confidential 
communication from a source node to its destination in the wireless cooperative communication network as 
studied in [9, 10]. In various system network such as cooperative relaying, jamming, multiuser scheduling, 
multiple-input-multiple-output (MIMO) with using PLS is considered in [11-15]. 

In this paper, we have proposed and investigated the Jammer Against Eavesdropper In Half-Duplex 
(HD) Energy Harvesting (EH) Cooperative Relaying Networks with one source node, one destination node, 
one intermediate relay, and in the presence of the Jammer Against Eavesdropper. We have analyzed the 
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system performance in terms of the integral form expression of secrecy outage probability. In addition, we 

have investigated the effect of source rate, time switching factor, energy coefficient, and the ratio P,/No on 

the system performance. Finally, all the mathematical, analytical expressions are verified by Monte Carlo 

simulation, and the analytical results match well with the simulation ones to convince the correctness of the 

analytical expression. The main contribution of this research can be formulated as the following 

a. The integral form of expression of secrecy outage probability is derived. 

b. The effect of source rate, time switching factor, energy coefficient, and the ratio Ps/NO on the system 
performance is demonstrated. 

c. The Monte Carlo simulation is conducted to verify the correctness of the analytical expressions. 


2. SYSTEM MODEL 

We consider a communication scenario with the help of a friendly jammer (J), as illustrated in 
Figure 1. The EH and IT processing are proposed in Figure 2. In this scheme, T is the block time in which the 
source fully transmits the information data to the destination. In the first interval time (aT), the R and J 
harvest energy from the S signal, where a is the time switching factor a € (0, 1) In the two remaining 
intervals time (1-a)T/2, the S and R node transfer information to the D node [5-7, 16-24]. 





°. 


S= a 





fsp hiro 
S eS 
— EH Phase — — IT Phase aT (1-æ)T/2 (1-æ)T/2 
Figure 1. System model Figure 2. The EH and IT phases. 


2.1. Energy harvesting phase 
In the first phase, the source will supply the energy for both jammer and relay nodes. Hence, the 
harvested energy at the jammer and relay can be given as, respectively 


2 
ER = nPsaT |hsp\? (2) 


The average transmitted power at the jammer and relay nodes can be obtained from (1) and (2), 
respectively 





_ Ey _ nPsat|hs | _ 2 
J (4-a)T/2— (1-@)T/2 kPs|hs,| (3) 
Pr = kPs|hsel? (4) 
Where k = 2 
1-a 


2.2. Information transmission phase 
In the second phase, the received signal at the relay can be rewritten as 


Vr = hspXs +N, (5) 


Where hsp is the channel gain of S-R link, x, is the transmitted signal from source and n, is additive 
white Gaussian noise (AWGN) with variance No and E{|x,|*} = Ps which E{e} is expectation operator. 
In the third phase, the received signal at the destination can be given by 


Ya =hppXy + Nna (6) 
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Where hrp is the channel gain of R-D link, x, is the transmitted signal from relay and ng is (AWGN) 
with variance No and E{|x,|*} = Pp 

Here, we consider amplify and forward (AF) mode at the relay. Hence, the amplify factor can be 
given as 


x PR 
an mR 7 
yr Ps|hspi2+No ( ) 


Substituting (7) into (6), we have: 


Ya = hrpôyr = hrpôlhsrXs +n] + na = hsrhrpôXs + hppdn, + na (8) 
signal noise 


From (8), the end to end signal to noise ratio (SNR) of S-R-D link can be calculated as 


y _ E{|signal|?} _ |hsrl|?lhrpl?8?Ps _ Ihsrl?lhRDl?Ps 
SRD 7 „o2 252 N 
E{|noise|*} lARD|“5“Not+No larDl? No +52 


(9) 


After doing some algebra and using the fact that No<<Pr, (9) can be rewritten as 


_ lnsri*lhrpl*PsPR 
YsrD = PrlhRDl? No+NoPslhsrl? (10) 
Substituting (3) and (4) into (10), the end to end SNR can be reformulated as 


_ k@lhsr|*lhepl? _ koxXY 
YSRD klhpp|2+1 kY+1 


(11) 


Where ® = A = teal E= Ml 


3. OUTAGE PROBABILITY (OP) 
The received signal at the eavesdropper can be given by 


Where hsg is the channel gain of S-E link and ng is AWGN with variance No and E flal} =P, 


To protect information from being eavesdropped by E, the friendly jammer J performs jamming. 
Hence, the SNR at the eavesdropper can be expressed as 


lnsgl*P 
Ye = OS = (13) 


Py\hyE| +No 


Substituting (3) and (4) into (13), we have: 





Plhspl? pZ 
Se a 14 
4 Kolhsj| |hje| +1 KOV+1 a 
3 2 2 
Where Z = [Asg] JV = [as] DA 
Next, the channel capacity of S-R-D and of eavesdropper can be obtained as, respectively 


d= 


Rsrp = TA log, 1+Ysrp) (15) 


1— 
Re = —" log, (1 + yg) (16) 


Lemma 1. The cumulative distribution function of V can be computed as 
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Utilizing the result in [25], the CDF of X and Y can be shown as the below (18) 


F(a) =1- 2 As Apak,(2,/AsAjea) 
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(17) 


(18) 


where K,,(¢)is the modified Bessel function of the second kind and v" order and Asz, Ajg are mean of random 


variables (RVs) hs Hee f respectively. 





From S the probability density function (PDF) of V can be calculated as, after applying the 





following formula 2 =n) = —K,_,(z) - “Ka (z) 


fv(@) = a) = 2./As Ape Ko(2,/As Aye a) 


3.1. Secrecy Outage probability (SOP) 
The Secrecy capacity of the system can be defined as 


R = max(0, Rsrp — Re) = tog: (e) 


Where [x]* = max(0, x) 
The SOP can be formulated by 


1+ 
SOP = Pr(R < R,) = Pr (ee 


a < p) = Pr(¥srp < P + PYgr — 1) 
= fo FysroCP + PVE — 1lYr = x) fy (dx 
Where p = 2*s and R, is source rate. 


In order to calculate the probability in (21), we have to find F, -pp (a) and f,,,(b) 
At first, we have: 


Fy spp (4) = Pr(Vsrp < a) 


Substituting (11) into (22), (22) can be rewritten as 





F A Pr(x <= +>) 
ne O = PE yaa j P KOY 


=| r (5+ -5y lY = y) Ody 


f l-fini eo 
=1-App exp (- terc) i exp (- ose) x exp(—Appy) dx 





Where Asp, App are the mean of RVs|/cp|*, |Agp|7, respectively. 
Apply equation [3.324,1] of the table of integral, we can obtain as followings 


Fysap(@) = 1 — 2 exp (—*882) 88488 x K, (2 ‘aun 


Next, combine with (14), the CDF of yp can be given by 
pZ 
F, (x) = Pr(yg < x) = Pr ( 


KPV +1 = 
oO xX 
= f Fz (E+ kxv|V = y) fray 
Applying the result from (19), equation (25) can be computed as 


x) = Pr(Z <Š + kav) 
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(20) 


(21) 


(22) 


(23) 


(24) 


(25) 
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F, (x) =? lAs,Aye [ ( — exp |-Asz (= + kxy)]} x Ko (2 Ass Arey) dy 


A oO 
= 1 — 24 AsjÀjg exp (- SE Í exp(—kxàsgy) x Ko(2/AsAjey)dy (26) 





Where Ac, is the mean of RV|hsgl? 
Applying Taylor series as follows 


(—kxAspy)" o (-1)"(kxAsp)” 
a l a E (27) 


exp(—kxAsey) = Xn=0 
Substituting (27) into (26), we have: 


A o (=1)”(kxàsg)” ræ 
Ep) = 1 — 2 fisije exp (— SZ) x Digg E [yt x Ko(2 As Ayey)dy(28) 


n! 


By changing variable t = Jy and then applying equation [6.561,16] of table of integral, equation 
(28) can be reformulated as 


4 yo (-1)"(kxAse)" __ ÀSEX 2 
F(x) =1 Daoa years X exp ( =) X [P(n + 1)] (29) 


From (30), the PDF of yç can be obtained as 





OF,(x) Wo (D (kx)" Asp)" t! Aspx 
f = =) x exp h-E) x Irn + OP 
* am MOAsAye)  ® 
_1)71 nz,n-1 
— Yip ae X exp (- erx) x [Fin + 1] (30) 


(n-1)!(AsJåjJeE) 


Applying results from (23) and (30) for (21), finally the SOP can be claimed by 


o) Asplpt+px-1] AspArRp|pt+px-1| AspARp|p+px-1]| 
SOP = f; fa — 2 exp (- ae) “setae PIPE x K, (2 jatar) f, (œ)dx 
(31) 


4. NUMERICAL RESULTS AND DISCUSSION 

In this section, the effect of source rate, time switching factor, energy coefficient, and the ratio P,/No 
on the system performance is investigated. The SOP versus the source rate Rs is plotted in Figure 3 with the 
main system parameters as d= 5 Db, n=0.8 and a=0.25, 0.5, 0.85. As shown in Figure 3, the SOP 
significantly increases while the source rate increases from 0 to 6 bps/Hz. Figure 4 presents the influence of o 
on the system SOP. In Figure 4, we set n=0.8, a=0.5 and Rs= 0.5, 1.0, 1.5 bps/Hz, respectively. From the 
results, we can see that the system SOP crucially falls down with the rising of o from 0 to 30 dB. In Figure 3 
and Figure 4, the simulation results match well with the analytical ones for verifying the correctness of the 
system performance analysis. 

The effect of the energy efficiency coefficient on the system SOP is illustrated in Figure 5, in which 
we set a=0.5, Rs=0.5 bps/Hz, o=1, 5,10 dB, respectively. Here, we can state that the system SOP has fallen 
down when n varies from 0 to 1. Furthermore, the system SOP versus time switching factor a is drawn in 
Figure 6. We set d=5 dB, Rs=0.5 bps/Hz, and n=0.25, 0.5, 1, respectively. The same as the above case, the 
system SOP decreases massively with the rising of a from 0 to 1. As shown in Figure 5 and Figure 6, the 
analytical and the simulation results are the same to convince the analytical system performance analysis. 
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SOP versus RA with ®=5dB,7=0.8 
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Figure 3. SOP versus Rs. 


SOP versus 77 with a=0.5,R.=0.5 bps/Hz 
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Figure 5. SOP versus n. 


5. CONCLUSION 
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SOP versus © with 7=0.8,a=0.5 
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Figure 4. SOP versus 9. 
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Figure 6. SOP versus a. 


In this paper, we have investigated the HD EH Cooperative Relaying Networks in the presence of 
the Jammer Against Eavesdropper. We have analyzed the system performance in terms of the integral form 
expression of secrecy outage probability. In addition, we have investigated the effect of source rate, time 
switching factor, energy coefficient, and the ratio Ps/No on the system performance. Finally, all the 
mathematical, analytical expressions are verified by Monte Carlo simulation, and the analytical results match 
well with the simulation ones to convince the correctness of the analytical expression. 
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